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ABSTRACT 

Do some Wolf-Rayet stars owe their strong winds to something else besides radiation 
pressure? The answer to this question is still not entirely obvious, especially in cer- 
tain Wolf-Rayet subclasses, mainly WN8 and WC9. Both of these types of Wolf-Rayet 
stars are thought to be highly variable, as suggested by observations, possibly due to 
pulsations. However, only the WN8 stars have so far been vigorously and systemat- 
ically investigated for variability. We present here the results of a systematic survey 
during 3 consecutive weeks of 19 Galactic WC9 stars and 1 WC8 star for photometric 
variability in two optical bands, V and I. Of particular interest are the correlated vari- 
ations in brightness and colour index in the context of carbon-dust formation, which 
occurs frequently in WC9 and some WC8 stars. In the most variable case, WR76, we 
used this information to derive a typical dust grain size of ~ 0.1 /im. However, most 
photometric variations occur at surprisingly low levels and in fact almost half of our 
sample shows no significant variability at all above the instrumental level (ct ~ 0.005 
- 0.01 mag). 

Key words: stars: Wolf-Rayet - circumstellar matter - stars: variable: other - stars: 
carbon - stars: individual (WR76) 



1 INTRODUCTION 

• Population I Wolf-Rayet (WR) stars are the descendants 
of the most massive stars on the main sequence, namely O 
stars. They present the signatures of strong, dense winds, 
in which typical terminal velocities are of order 10'' km s~^. 
WR spectra come in two sequences, those with strong lines of 
helium and nitrogen (WN subtypes) and those with strong 
lines of helium, carbon and oxygen (WC and WO). WC9 
are the coolest stars among population I WC subtypes, with 
spectroscopic temperatures of order 50000 K at the stellar 
(hydrostatic) radi us, deduced from non-LTE, clumped, line- 
blanketed models ijDessart et al" I2OOOI : ICrowther et aLll2002l : 
iBarniske et al.ll2006[ ). The particularity of WC9 stars is their 
frequent very bright lumin osity in the infrared t hat is associ- 
ated with dust formation ( Willi ams et al.l|l987l ). Some WC9 
stars al so show strong s pectroscopic and photometric vari- 
ability llCrowtheij|l997l : lEenens fc Corrall bood : iKato et al.l 



by recent MOST satellite data l|Moffat et al.l 120081 ). A few 
other WC9 stars have been monitored in the optical by 
IVeenI l|2000l ). revealing behaviour that reminds one of the 
all-yariable WN8 stars, which typically shows a level of 
photom etric variability of g ~ 0.016 to 0.05 mag i n the 
visible llLamontame fc Moffaj Il987l: lAntokhin et al1ll995l : 



iMarchenko et al.Ml998l : iLefevre et al.l |2005^ ■ Photometric 
variability of other WC9 stars both in the infrared and in 
the visible can be explained by episodic or periodic dust for- 
mation (NIR emission and optical absorption) associated in 
most cases with binarity (Veen et al. 19981: Williams Il996l : 



l2002l : lLepinelll996l ). The brightest WC9 star, WR103, was 
observed to show periodic behaviour , suggesting a poss l 



IWilUams fc van der Huchtll2000l : IWilliams et allboosl ). 

The goal of this study is to provide a systematic, uni- 
form survey of WC9-stars using precision optical photome- 
try in an attempt to answer the questions raised by previous 
studies, namely: Do all WC9 stars show variability like that 
of WN8 stars? and Is variability associated with dust for- 
mation or with pulsations or both? 



ble binary (jlsserstedt fc MoffatI Il98ll: iMoiflat et al 



although later more intense data (|Veen et all Il999l ) sue- 2 OBSERVATIONS 



gest that stellar pulsations are at play, as also confirmed 
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We selected a complete sample of 19 WC9 stars and 1 WC8 
star in the southern sky from the "Vllt h Catalogue of Galac- 
tic W olf-Rayet Stars" and its annex (|van der Huchtll200ll . 
I2OO6I ) with V-band magnitude from 12 to 16, in order not to 
go below the minimum practical exposure time of the CCD 
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and to stay below a reasonable maximum exposure time of 
600 seconds. Observations were made using the Swope 1- 
meter telescope at Las Campanas (Chile) between UT 2007 
June 8 and July 7. We used the direct CCD camera SITe#3 
and the two optical bands V and I in order to be sensitive 
to both magnitude and colour variations, which are of par- 
ticular interest to constrain the local dust properties. Our 
field of view is 8.5 x 10 arcminutes with a pixel size of 0.435" . 
We chose net exposure times such that we obtained a good 
signal-to-noise ratio (at least ~300) per observation with one 
or two observations per clear night for each star. Standard 
star fields were not observed, since we are only interested in 
differential photometry. 



3 DATA REDUCTION 

Extraction of stellar magnitudes from the images was made 
in the usual way, using nightly flat-field images in each 
band, nightly bias images and overscan pixels in each im- 
age. We then chose to perform aperture photometry on our 
images, which is more efficient and precise than the alter- 
native point-spread extraction for well-isolated stars in the 
field. We first interpolated all the images to match a refer- 
ence frame so that the target star always has the same pixel 
osition. Th is was done with the ISIS routine interp.csh 
Alardll2000l ). We then chose an adequate number of refer- 
ence stars (Niisf) for which we measured the instrumental 
flux F. To do so, we used the IDL function APER with 
a fixed aperture radius R — 6 pixels (2.610") in I band 
and R = 8 pixels (3.480") in V band, which correspond to 
2xFWHM for typical images, together with a sky annulus 
of radii [1.5 x R,2.5 x R] . The instrumental flux is then 
computed as 



sky 



where P(r) is the pixel value at radius r and Fsky is the 
mean fiux computed in the sky annulus. The differential 
magnitude AMi of each star in the field is finally given by : 



AM, 



-2.5i 



F, 



We then sorted the stars in order of variance and chose the 
Nn^f least variable among them to be the new reference 
stars. We then re-performed the process until it converged. 
Nnsf is chosen in order to minimize the variance of the ref- 
erence stars and ranges from 10 for a very diffuse field to 
50 or more for a crowded field. The photometric precision 
reached goes from (j(AM) = 0.005 to 0.012 magnitude de- 
pending on the quality of the images, essentially determined 
by the average airmass at which the star could be observed. 



4 RESULTS 

We found 12 WC9 stars to be variable when compared to ref- 
erence stars of similar magnitude : WR76, WR48b, WR59, 
WR75C, WR81, WR88, WR121, WR77t, WR80, WR92, 
WR96 and WR119 (FigH] and [H. Their variability ranges 
from a = 0.012 to 0.06 mag, the last value referring to the 
by far most variable star of the sample, WR76. Aside from 
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Figure 5. V band differential magnitude versus V-I differential 
colour index for WR76, by far the most variable star of our survey. 
The expected curves for Rayleigh scattering (valid for dust grains 
smaller than A^/IO ~ 0.05 /xm) and grey scattering (valid for 
dust grains larger than Xy — 0.5 /.tm) are shown for comparison. 
The Rayleigh slope has been computed for Mie scattering with 
Qext(X,a) oc A""* (see text). 



this star, the level of variability of our sample is quite small 
(going from a 0.012 to 0.02 mag) rela tive to the "violent" 
WR103 described bv lVeen et al.l (|1999l ). who confirm a value 
of (j(AM) of ~0.035 mag in V band (0.05 mag p eak to val- 
ley) initially found by llsserstedt fc MoffatI ( *1981) and con- 
firmed bv iMoffat et all (119861 ) and lBalonaet al . (1989). The 
timescale of the (stochastic) variability of WR76 lies in the 
range of ~ a day to several days, like that seen in WR103. 

Eight stars were not detected as variables above the 
instrumental value (FigO and |4]). Among these stars we 
estimated the upper bound for ay to be around 0.01 for 
WR106, WR95 and WR65 and around 0.005 for WR53, 
WR73, WR75a, WR75b and WR117. 



5 DUST FORMATION 

We searched for a correlation of V vs. (V-I) but only found 
one for WR76, which is by far the most variable star of 
our survey (FigO. The computed correlation coefficient for 
this star is R=0.92 and the Spearman rank index, which 
expresses the probability that the two variables are un- 
correlated, is 2.68 x 10"* . We observe that WR76 be- 
comes redder when the V extinction increases (the star 
gets fainter). This is consistent with circumstellar extinc- 
tion by small dust grains, which is more efficient in V band 
than in I band. Moreover, WR76 is known to have a very 
high dust formation rate measured from infrared photome- 
try l| Williams et al.lll987h : its dust formation rate is by far 
the highest among our sample. Nevertheless, the correlation 
between the dust formation rate and the level of variabil- 
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Figure 1. Dispersion plots for the fields containing the 12 variable WC9 stars: WR76, WR48b, WR59, WR75c, WR81, WR88, WR121, 
WR77t, WR80, WR92, WR96 and WR119. These are for the V band except for WR77t, which has been observed only in I band. The 
reference stars are indicated by triangles and the Wolf-Rayet by a square. On the x-axis is the instrumental magnitude of the stars 
measured with arbitrary zero point from an image chosen as reference. On the y-axis is the standard deviation in magnitudes of AM 
with time (cr). 



ity of the star is completely dominated by this star, as seen 
in Fig|6] For the remaining stars, the correlation is quite 
poor. In particular, we found variable stars of low amplitude 
with no dust formation detected (WR75c, WR81, WR88, 
WR92) together with non variable stars known to produce 
dust because of their strong infrared excess (WR53, WR73, 
WR117). This might be at least partly due to the shorter 
timescale of our variability search compared to that for typ- 
ical dust ejection. 

As some non-dusty WC9 stars show low-level variabil- 
ity, presumably some of that observed from WR76 may also 
be intrinsic. However, on account of its high dust formation 
rate, variable extinction is thought to play a dominant role. 

Estimations of the dust formation rates were not avail- 
able for all our target stars, although the analysis of this 
correlation does show quite clearly that the production of 
dust cannot explain the behaviour of all WC9 stars. Like in 
the strongly varia ble WC9 star WR103 (|Veen et al.lll999l : 
I Moffat et al.llioO^ '). we believe that pulsations are likely to 



play an important role in some cases, sometimes combined 
with the effects of dust formation. 



6 DUST PROPERTIES 

It is known that dust in WC9 stars is essentially composed 
of amorp hous carbon grain s (|Williams et al.l Il987l : IZubkol 
1998; Ma rchenko et al.ll2q02f ). We therefore used Mie theory 
(jBohren fc Huffman|[l98^ ) and a model of monosize carbon 
spheres, with optical constants from lZubko et al. (1996), in 
order to interpret the measured slope of Fig[S] as a domi- 
nant size in the dust size distribution. To do so, we calcu- 
lated the weighted extinction coefficients Qext{^,a) (func- 
tion of the wavelength A and the dust sphere radius a) as in 
iMarchenko et al.ll|2003D . the weights being given by W{X) = 
T{\)F{\), where T(A) is the transparency of the correspond- 
ing filter (V or I) and F(\) is the energy distribution of 
WR7 6 in the optical and the near infrared (jWilliams et al.l 
Il987t ). We then used the fact that the extinction, A{\), is 
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Figure 2. V and I lightcurves for the 12 variable WC9 stars compared with a reference star of approximately the same ma gnitude in V. a 
is the rms st andard deviation of AMag in time. When available, Vb and Iq are the apparent magnitudes of the stars from lvan der Huchtl 
||2001| . 120061 ') . When not available, Vq and Iq are estimated values deduced from an extrapolation of what was available in the catalogue. 
MJD is for modified Julian day. 



directly proportional to Qext- Finally, A{X) is linked to the 
variations of magnitude that we interpret as due to an in- 
crease of the extinction by the dust {A{X) oc AMag ), The 
result of our model is shown Fig[7]for a ranging from 0,0001 
to 0.2 nm. Our estimated dust size, corresponding to a slope 
of 4.47±0,60, is a = 0.126±0,003 fim. This value of a is quite 
close to the estimation of dust size in WR140 (WC7pd-F05) 
during its periastron t ransit in 2001: a — O.OBQIq ooi M"^ 
l|Marchenko et al.l [20031 ') but is somewhat lower than other 
estimations in W R112 (WC9d-FOB?) aiid 'WR118 (WC9d) : 
g = 0.4 - 1.0 urn llChiar fc TielensllioOll : lYudin et al.ll200ll : 
iMarchenko et al.ll2002l ) . 



7 CONCLUSIONS 

The most surprising result of this study is the large fraction 
of WC8/9 stars that are not variable above the a ~ 0.01 
mag level. Even among the 12 that do vary, it is mostly 
at a relatively modest level, cr < 0.02 mag, in contrast to 
their highly variable WN cousins, the WN8 stars. We see 



a significant variation in both V and I filters for only one 
object, 'WR76, which we ascribe to variable circumstellar 
dust extinction. 'We estimated the dominant dust grain size 
to be close to 0.1 fj,m, which is ~ an order of magnitude 
smaller than estimated in two other WC9d stars. 
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Figure 3. Dispersion plots for the fields containing the 7 non-variable WC9 stars WR106, WR95, WR65, WR73, WR75a, WR75b and 
WR117, and the non-variable WC8 star WR53 as in Fig.4. These are all for the V band. 
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Figure 4. V and I lightcurves for the 8 non-variable WC9 stars compared with a reference star of approximately the same magnitude 
in V . 
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Star name Spectral Type av of WR ay of typical reference star aj of WR crj of typical reference star Slope 



WR76 


WC9d 


0.060 


0.004 


0.051 


0.006 


4.47±0.60 


WR48b 


WC9d 


0.018 


0.005 


0.017 


0.005 


- 


WR59 


WC9d 


0.012 


0.005 


0.017 


0.007 


- 


WR75C 


WC9 


0.015 


0.005 


0.017 


0.005 


- 


WR81 


WC9 


0.016 


0.010 


0.021 


0.010 


- 


WR88 


WC9 


0.017 


0.005 


0.021 


0.006 


- 


WR121 


WC9d 


0.019 


0.008 


0.020 


0.009 


- 


WR77t 


WC9d 


- 


- 


0.017 


0.005 


- 


WR80 


WC9d 


0.012 


0.006 


0.011 


0.007 




WR92 


WC9 


0.013 


0.008 


0.014 


0.007 




WR96 


WC9d 


0.013 


0.004 


0.011 


0.007 




WR119 


WC9d 


0.014 


0.006 


0.012 


0.005 




WR106 


WC9d 


0.012 


0.010 


0.012 


0.008 




WR95 


WC9d 


0.009 


0.007 


0.012 


0.005 




WR65 


WC9 


0.001 


0.004 


0.010 


0.007 




WR75a 


WC9 


0.004 


0.005 


0.007 


0.007 




WR53 


WC8d 


0.005 


0.007 


0.007 


0.006 




WR73 


WC9d 


0.005 


0.005 


0.005 


0.005 




WR117 


WC9d 


0.005 


0.007 


0.006 


0.008 




WR75b 


WC9 


0.008 


0.005 


0.009 


0.005 





Table 1. Measured V and I variability compared to typical values for the reference stars. The latter is estimated from the rms mean of 
the reference stars a. 'Slope' is the slope of the (V) vs (V-I) plot (Fig. 1), which was only significant for the most variable star of the 
survey, WR76. 




Figure 6. Our measured level of photometric variability a[V) 
versus the estimated dust formation rate by IWilliams et al.l 
lll987l '). We show the range of values for the instrumental level, 
which is determined, for each field, by the cr of the reference star 
of closest magnitude to the WC9 star. The dust formation rate 
estimates are based on infrared photometry. We see that any cor- 
relation between the level of variability and the dust formation 
rate of the WC9 stars is completely dominated by only one star, 
WR76, combined with a cloud of points at low variability level. 
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Figure 7. Model prediction of the V vs. (V-I) slope for WR76. 
We used Mie theory for mono-size d spherical grains p lus optical 
constants for amorphous carbon bv lZubko et al.l lll996l) . The best 
model to reproduce our slope of 4.47 ± 0.60 is a dust size of a = 
0.126 ± 0.003 ^lm 
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